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We present a framework for the quantum enhanced estimation of multiple parameters corresponding to non¬ 
commuting unitary generators. Our formalism provides a recipe for the simultaneous estimation of all three 
components of a magnetic field. We propose a probe state that surpasses the precision of estimating the three 
components individually and discuss measurements that come close to attaining the quantum limit. Our study 
also reveals that too much quantum entanglement may be detrimental to attaining the Heisenberg scaling in 
quantum metrology. 


Introduction: As the elementary theory of nature, quantum 
mechanics sets the fundamental limit to precision of param¬ 
eter estimation. On the flip side, quantum resources enable 
the estimation of parameters with a precision surpassing that 
set by classical physics. This is the basis of the field of quan¬ 
tum enhanced sensing and metrology, and has been studied in 
great depth both theoretically and experimentally [1-4]. Al¬ 
though most of these investigations have largely focussed on 
the estimation of a single phase parameter, more recently at¬ 
tention has been cast on the quantum-enhanced estimation of 
multiple parameters simultaneously [5-13], and some early 
experiments have already been performed [14]. 

The motivations for studying quantum-enhanced multi¬ 
parameter estimation are manifold: First, while single phase 
estimation captures a wide range of scenarios [15], high-level 
applications such as microscopy, optical, electromagnetic, or 
gravitational field imaging, and spectroscopy intrinsically in¬ 
volve multiple parameters that should be estimated simultane¬ 
ously. Secondly, while the quantum-enhanced limit for indi¬ 
vidual phase estimation can always be attained [16, 17], the 
measurements required to attain the quantum-enhanced limit 
for multiple parameters need not necessarily commute. This 
makes multi-parameter quantum-enhanced sensing a very in¬ 
teresting scenario for studying the limits of quantum mea¬ 
surements [6, 7]. Finally, multi-parameter quantum-enhanced 
sensing provides a novel paradigm for investigating the infor¬ 
mation processing capabilities of multi-partite or multi-mode 
quantum correlated states and measurements. 

In this work, we study the problem of estimating a multi¬ 
dimensional field using a fixed number of particles. We first 
show that for a uniform field, the quantum enhancement to the 
precision of estimation is provided entirely by the two-particle 
reduced density matrix of the system, and that the attainabil¬ 
ity of the quantum enhancement is solely determined by the 
one-body reductions of the probe state. We apply our meth¬ 
ods to the simultaneous estimation of all the components of a 
magnetic field in three dimensions and show that this can be 
about three times better than estimating the components indi¬ 
vidually [18-21]. Finally, we present a multi-partite quantum 
state achieving this advantage, and show how realistic mea¬ 
surements perform in attaining the multi-parameter quantum 
limit using matrix product state techniques [22-24]. 

Framework: We consider the estimation of parameters gov¬ 


erned by the Hamiltonian 

d 

Hiv) = E (1) 

k=l 

The parameters pk C R, k = 1,..., d, to be estimated are 
the coefficients of a set of (not necessarily commuting) gen¬ 
erators Hk . We assume that the Hk themselves do not depend 
on cp. In addition to estimating a field in multiple dimensions 
simultaneously in free space, materials or biological samples, 
this problem is equivalent to quantum-enhanced Hamiltonian 
tomography as it allows us to estimate unknown coefficients 
of the Hamiltonian in a suitable operator decomposition [25]. 
We note that earlier works have studied the estimation of pa¬ 
rameters corresponding to unitary channels from information 
geometry [26-28] and representation theory [29, 30] perspec¬ 
tives, and have shown a Heisenberg scaling in their estimation. 
A pure W-particle probe state \fj) acquires the parameters via 
the unitary transformation U{p>) = and we seek 

the best quantum strategy for the estimation of the param¬ 
eters from the evolved probe state = U{cp)\fj). The 
performance of an estimator of (p is quantified in terms of 
the covariance matrix Cov [(/::>]. The quantum Cramer-Rao 
bound [16, 17] is a lower bound to the covariance matrix in 
terms of the quantum Fisher information matrix (QFIM), thus 
yielding an ultimate limit on the best possible precision of any 
(unbiased) estimator. For every specific set of positive opera¬ 
tor valued measurements (POVM) {H^}, one finds [17] 

MCoy[^] > {n,})-i > X((^)-^ (2) 

where the first inequality is the classical and the second in¬ 
equality the quantum Cramer-Rao bound, respectively. Here, 
M is the number of times the overall experiment is repeated 

and Tk,i{(pytli}) = Tjn^VkP{n\(p)d^iP{n\‘^)/p{n\(p), 

k,l = 1,..., d, denotes the Fisher information matrix (FIM) 
determined by the probabilities p(n1 9 :?) = Fur¬ 

ther, Xk^i{(p) = Re[{fj(p\LkLi\'ipy:,)] is the QFIM, where, for 
pure probe states, the symmetric logarithmic derivative (SLD) 
Lk with respect to the parameter (/ 9 /c is determined hy Lk = 
2 [\9cpki^cp){^cp \ + \i^cp){dcpdPcpW for ah k = 1,.. . [17]. 

While the classical Cramer-Rao bound can always be satu¬ 
rated by, e.g., a maximum likelihood estimator [31], the quan¬ 
tum limit (i.e., the second inequality in Eqn. (2)) may not be 
attainable in general. In a single parameter setting, the optimal 
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measurements saturating the quantum Cramer-Rao bound are 
given by the projectors onto the eigenvectors of the SLD. In 
the multi-parameter setting, however, the SLDs may not com¬ 
mute in general, leading to tradeoffs for the precisions of the 
individual estimators [6, 7] . 

Formalism: For unitary time evolutions under Hamiltonians 
of the form of Eqn. (1) the QFIM can be expressed as the 
correlation matrix of the Hermitian operators defined by [32] 
(see Appendix A) 

Jo 

leading to (suppressing the parameter cp in the arguments 
henceforth) 

Ik,I = 4Re[(V’|ifci(|V’) - (4) 


We now restrict to the situation where the N particles evolve 
under the one-particle Hamiltonian = Ylk=i for 

n = 1,..., (where the h^, ^ are bounded), leading to the 
global Hamiltonian 

N d N d 

= = (5) 

n=l /c=l n=l /c=l 

With this, the operators defined in Eqn. (3) simplify to 
Akiv) = ELi 4"'> where 

Jo 

are Hermitian operators acting only on particle n. Incorporat¬ 
ing this in Eqn. (4) simplifies the QFIM to 


N 


Ik,i=^'Y^ Re - (V’l4"'lV’)(^|a|™'l^) 


n,m=l 


(7) 


= 4^Re[Tr[pWaNa{”l] -Tr[^WaM] Tr[pWaW]] +4^ Re[Tr0 a{“'] -Tr[pWaM] Tr[^HaH] 

n n^m 


where = Tr\|^ [j?/;) (t/^I] denotes the reduced density 

matrix to sub-systems n, m. The first sum has N terms, and 
the second has 0{N‘^) terms. The latter points to the origin of 
the quadratic scaling in N in quantum-enhanced metrology, 
and shows that quantum advantage in metrology can be no 
more than quadratic in the considered scenario of Eqn. (5) if 
the variance of the estimator is the quantifier of precision. 

For estimating a uniform field, the phase parameters are 
identical across the system (although corresponding to non¬ 
commuting generators) and we can restrict ourselves to per- 
mutationally invariant quantum states, i.e., states that are in¬ 
variant under any exchange of its constituents. The one- 
and two-particle reduced density matrices are then given by 
gH = ^[1] and for all n, m, respectively. Under 

the restriction of permutationally invariant states, the QFIM 
simplifies to 

I = 4NI^^^ + 4N{N - 1)2:[2] ^ (8) 

where 

4'i = Re[Tr[pWafca;]] - Tr[pWafe] Tr [^Wa^] (9) 

only depends on the one-particle reduced density matrix and 

xg = Tr[^p]^^ ^ _ Tr^ilafc] Tr[pWa;] (10) 

depends on the two-particle reduced density matrix. 

Eqn. (8) highlights several interesting physical aspects of 
quantum-enhanced metrology: First, note that can be 
bounded independently of . This immediately shows that 


the archetypal quadratic scaling of quantum-enhanced sens¬ 
ing arises solely from the two-particle reduced terms. For 
instance, let the probe state he \ 2 p) = 10)®^, i.e., permu¬ 
tationally invariant and separable. Then, (g) 

such that = 0 and the QFIM only scales linearly in N, 
i.e., X = Thus, Eqn. (8) implies that in permutation- 

ally invariant systems quantum correlations are necessary for 
achieving a quadratic scaling in the number of probe states 
N —the so-called Heisenberg scaling. Note that the latter rea¬ 
soning also applies to quantum states that are not permuta¬ 
tionally invariant, as can be seen by Eqn. (7). Further, for 
probe states of the form \^jj) = the QFIM satisfies 

rank[X] < 2{D — 1) where D is the dimension of the lo¬ 
cal Hilbert space (e.g., D = 2 for two-level systems, see 
Appendix B for details) such that if the number of param¬ 
eters exceeds 2{D — 1), i.e., d > 2(7) — 1), a simultane¬ 
ous estimation of all parameters necessarily fails due to a 
lack of information for all parameters in the QFIM. Finally, 
if both the one- and two-particle reduced states are maximally 
mixed, the Heisenberg scaling is lost. To see this, note that 
gi^] = I 2/2 (where 1/c is the k x k identity matrix) implies 
= Tr (g) di] — Tr [di] Tr [a/]/4, which vanishes 

if ^[2] = I 4 / 4 . This is an example where too much entangle¬ 
ment harms the quantum advantage of exploiting N particles 
in parallel. 

Attaining the quantum limit: Saturating the quantum Cramer- 
Rao bound and attaining the QFIM is the next important part 
of quantum-enhanced sensing. This is particularly interesting 
for multi-parameter estimation since the SLDs corresponding 
to the different parameters need not commute. We show in 
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Appendix C that for a purely unitary evolution, the QFIM is 
saturated if the expectation value of the commutator of the 
SLDs vanishes for all pairs [28], i.e.. 


(I 2 (8) I 2 + ^ <^/c)/4, respectively. Now, for the simultane¬ 

ous estimation of all three components, an obvious candidate 
is 


{tp^lLkLi - LiLk\^^) = SilmliiplAkAilip)] = 0. (11) 

For permutation invariant systems, this reduces to 
8iA^Im[Tr [p[^]d/cd/]] = 0 for all k,l. It is interesting 
to note that while the quantum-enhanced scaling is governed 
entirely by the two-particle reduced density matrices (see 
Eqn. (7)), the attainability of this bound is determined 
solely by the one-particle term (for a general proof, see 
Appendix C). The expectation value vanishes, for instance, 
for permutationally invariant pure probe states with 
^[1] = I 2/2 [33]. This is a sufficient but not necessary 
condition for the expectation of the commutator to vanish and 
gives a rather simple mathematical condition for the quantum 
Cramer-Rao bound to be saturated. It is an instance of the 
local suppression of the non-commutativity of the generators 
using quantum correlations [26] . 

More generally, when the expectation values of all commu¬ 
tators of the SLDs vanish and the QFIM is of full rank, the 
eigenvectors of the d distinct SLDs lie in a subspace of dimen¬ 
sion d + 1 allowing for the construction of a POVM that satu¬ 
rates the quantum Cramer-Rao bound. We prove this assertion 
in Appendix C and, further, provide a procedure for construct¬ 
ing such a POVM that saturates the quantum Cramer-Rao 
bound. Note that for commuting generators { 2 p\AkAi\ 2 p) G R, 
such that the quantum Cramer-Rao bound can always be sat¬ 
urated given the QFIM is not rank deficient (see also [28]). 
Estimating a magnetic field in three dimensions: We now ap¬ 
ply our formalism to the task of estimating the components 
of a magnetic field in three dimensions simultaneously using 
two-level systems. Potential systems could include trapped 
ions, NV centres, or doped spins in semiconductors [34-38]. 
The Hamilton operator for this system is given by (with d = 3, 
see Appendix D for a discussion of d > 3) 

33 3 

h = fi - B = ^ fikBk = y] := ipkak, (12) 

/c=l k=l k=l 

where the magnetic moment fik = p^dk /2 is proportional to 
the spin, {dk} denotes the unnormalised Pauli operators, and 
(pk = l^Bk/2. To develop the intuition for estimating the 
magnetic field in three dimensions simultaneously, we start 
with the estimation of a magnetic field pointing solely along 
one of the specific directions X, Y, or Z. It is well known that 
a GHZ-type state (see Appendix E) 

\^k) = {\^t)'^^+ \^k)^^)/V2 (13) 

achieves the quantum Cramer-Rao bound, where 1) is the 
eigenvector of the Pauli operator dk corresponding to the 
eigenvalue ±1 {k = 1,2,3 corresponding to the X, Y, and 
Z direction). These states are permutationally invariant with 
one- and two-particle reduced density matrices = I 2/2 
and = (10^, 0^1 + 10^ , |)/2 = 


\^) = 1 ^ 1 ) + 1 ^ 2 ) + l^s)) , (14) 

where M is the normalisation constant and {(5/c} adjustable 
local phases. Now, for even N and appropriate 5^ the quan¬ 
tum Cramer-Rao bound is achievable with this state since 
^[1] = 12/2 [39]. It remains to show that it is also capable 
of achieving the Heisenberg scaling. To simplify our calcula¬ 
tions, we henceforth restrict ourselves to systems that consist 
of = 8n, n G N, particles (and Sk = 0 for all k) but note 
that this is no limitation of our model as indicated by the nu¬ 
merical simulations presented below. For N = 8n, n G N, the 
two-particle reduced density matrix of 1 7 /^) is an equal mixture 
of the GHZ-type states in all directions and given by 

= 3 = 4^2 ® 12 + y^CTfc(gxTfe. (15) 
^k=l k=l 

For any other N, we show in Appendix F that the difference 
from the form of in Eqn. (15) is exponentially small in N. 
For a probe state with marginals = I 2/2 and given 
above, the QFIM is (see Appendix G and Ref. [27] which 
shows the same scaling) 

Bk,i = ^N{N + 2) [{l-smc^[^])rikm+h,I sinc^l^]] , (16) 

where sinc[^] = sin[^]/^ with ^ = ^/Pl + P 2 -^ ^3 and 
Vk = + ^2 + ^3 Note that, in the limit 

of pk ^0 for k = 1,2,3, the QFIM is diagonal, i.e., 
Ik,I = ^N{N 2)Sk,I. As our quantum Cramer-Rao bound 
can be saturated, the minimal total variance for estimating 
the three components of the magnetic field simultaneously is 
given by 


= Tr[Cov((^)] = Tr[X-i((^)]. 

k=l 


(17) 


Since the QFIM in Eqn. (16) is the sum of a rank one matrix 
and a rescaled identity, its eigenvalues can be read off directly 
as Ai = 4N{N + 2)/3 and Az.a = 4Ar(Ar + 2)sinc^K]/3 
leading to 


I A »im |2 ^ 3 + 6/sinc^[e] 
I I 4N{N + 2) 


(18) 


Let us now compare three different scenarios for the estima¬ 
tion of p: (i) A classical strategy of using only pure prod¬ 
uct states, (ii) a quantum strategy where the parameters are 
estimated individually (as depicted in Fig. 1 (a)), and (iii) 
the simultaneous estimation of the parameters, as shown in 
Fig. 1 (b). To obtain a fair comparison amongst (i)-(iii), we 
use exactly N particles to estimate all three cases. As for all 
scenarios the QFI depends on the values of pk, we restrict 
ourselves to a setting with pk 0, k = 1,...,3. For sce¬ 
nario (i), the strategy is to divide the set of N particles into 
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three blocks of length n = N/ 3 and, on the /c* block, to pre¬ 
pare a product state that allows for the estimation of ifk. This 
is due to the impossibility of estimating 3 parameters simul¬ 
taneously using a pure and permutationally invariant product 
state, as shown by the singularity of the QFIM (Appendix B 
shows that its rank is 2). The maximal QFI for each block 
(see Appendix E) is equal to Xk = n(Aniax(a/c) - Aniin(d/c))^ 
where Xmax/mm{^k) denotes the maximal/minimal eigenvalue 
of d/c such that 

(Amax(4) - An,in(aA:))^ = 4 [(1 - sinc^[^])7?fc + sinc^[^]] 

with (Amax(afe) - Ami„(afc))2 -> 4 for (yjfe 0, A: = 1,2,3. 
Further, A(/9^ = and thus we find for the individual esti¬ 
mation of all parameters using separable states 

= (19) 

k=l 

Secondly, for a quantum strategy exploiting entangled states 
where we estimate the parameters individually we, again, di¬ 
vide the chain of N particles into three blocks. Next, on the 
block, one prepares a GHZ-type state in the ak basis. Re¬ 
call that for each block X/. = ^^(Aniax(^/c) ~ Aniin(d/c))^ (see 
Appendix E) such that with Aip‘1 = l/X/^ one finds 

= Jj. (20) 

k=l 

Thirdly, for the simultaneous estimation of the parameters we 
have (see Eqn. (18)) 

Hence, it is possible to design quantum probes (those with 
a completely mixed one-particle reduced marginal and two- 
particle reduced marginal given by Eqn. (15)) for magnetic 
field estimation such that estimating the three components si¬ 
multaneously is about three times better than estimating them 
individually. Overall, < jAc^g^pp for 

all N > 3. This is illustrated in Eig. 1 (c) where the re¬ 
sults are obtained numerically using matrix product state tech¬ 
niques [22-24] to also account for system sizes N ^ Sn. Note 
that the 3-fold improvement is not proven to be optimal, al¬ 
though we expect it to be so. This is a more general behaviour 
than the findings of [8], where the generators were commuting 
while here they are not. 

Classical Fisher information: We have already shown (Ap¬ 
pendix C) that there is a POVM that achieves the multi¬ 
parameter quantum Cramer-Rao bound. The so-constructed 
POVM contains as one element the projector onto the time- 
evolved probe state, i.e., U (v^) | . While this set theoretically 
achieves the bound, it may not be very appealing from an ex¬ 
perimental perspective. Hence, let us finally discuss some re¬ 
alistic measurements. In particular, we consider two sets of 
POVMs: /c = 1,. .., 4, contains the three projectors 

= |^fc)(^,| with l^fe) = (1,/.+)®^ + \ci>-r^)/V2 


(a) 

individual 

estimation: 




(b) 

simultaneous 

estimation: 






Figure 1. (a): Individual estimation of the phases can be achieved by 
dividing the chain into sub-blocks, each of which is used to estimate 
only a single parameter. This is contrasted by the simultaneous esti¬ 
mation of the parameters depicted in (b). (c): Log-Log plot for the 
estimation of the three directions of a magnetic field with parameters 
<^i = 10“^,(y92 = 2x10“^, and <^3 = SxlQ-^. We show the QFIM 
for the three different scenarios described in the main text, as well as 
the result obtained for the FIM for the two considered POVMs. Note 
that for the QFIM results we computed the total variance for all V, 
while for the FIM results only for the values of N emphasised with 
a marker. 

together with the element guaranteeing normalisation H^^^ = 

1 — X]/e=i A Note that for even N and appropriate 5^ 
these operators indeed form a valid set of POVMs [40] . Fur- 
ther, H^ k = 1,..., 3, is determined solely by expectation 
values of simple Pauli strings, i.e., 

na = (l±ar)/6. 

Again, we use matrix product state techniques to compute the 
classical Fisher information for these POVMs. This enables us 
to analyse the scaling up to very large system sizes, as shown 
in Fig. 1 (c). Note that the precision for both these POVMs is 
best for N = Sn and the considered probe state. 

Conclusions: We have obtained the quantum limits for the 
simultaneous estimation of parameters corresponding to non¬ 
commuting unitary generators. We applied our methods to the 
simultaneous estimation of all three components of a mag¬ 
netic field in space. The results suggest that estimating the 
phases simultaneously improves the sensitivity by a factor 
of d = 3, in consonance with earlier results with commut¬ 
ing generators [8]. Future extensions of our results could in¬ 
clude, amongst others, a combination of commuting and non¬ 
commuting generators, and the inclusion of decoherence. An¬ 
other direction could be the search for optimal probe states 
and more tractable measurements for specific physical sys- 
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terns such as trapped ions or vacancy centres in diamond. Acknowledgements: This work was supported by the UK EP- 
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Appendix A: Unitary multi-parameter estimation: 
SLDs and QFIM 


In the first section of the appendix, we set out to find an ex¬ 
pression for the symmetric logarithmic derivatives (SLDs) to¬ 
gether with the quantum Eisher information matrix (QEIM) 
for the setting discussed in the main text. Eor this, we restrict 
ourselves to unitary channels where the to-be-estimated pa¬ 
rameters pk C R, k = 1 ,..., d, are the coefficients of a set of 
(not necessarily commuting) generators Hk, i.e., we consider 
unitaries of the form 

U{(f) = = e-‘ ^^=1 , (Al) 

where hI = Hkioiallk = 1 ,... ,d and <p G with = 
Pk. Eurther, note that the Hk do not depend on the parameters 
cp. Eor a pure probe state \ 2 p) and purely unitary evolution, the 
SLDs are given by [17] 


Lk — \ (A2) 

where {(p)]\'ip) denotes the partial deriva¬ 

tive of I'lpip) with respect to the parameter pk. Now, recall 
that [32] (see [18] for another application in quantum metrol¬ 
ogy) 


dpk 


■ = —1 


fdae ^ 

Jo 




dH[^) ,-i 
dpk 


iaH (if) 


(A3) 
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i.e., 

(A4) 

with the skew-Hermitian operator 

dfeM = -iife((^) = -i (A5) 

Jo 

and where we defined Ak{(p) = iOk{^)- With Eqns. (A2) 
and (A4) one finds 

L,=2?7[6fe|V’)(^| + |V’)(V’|Oi]f7^ 

- - - . (Ao) 

= 2iU[\^P){^\,Ak]U^, 

where [X,Y] denotes the commutator of the operators X and 
Y, respectively. Next, let us consider the QFIM. For unitary 
time evolutions it is given by [ 8 , 17] 

Ik,I (<P) = 4 - (d^^4’<pl4’<p) ■ 

With this, Eqn. (A4) allows us to write the QFIM in terms of 
the correlation matrix of the operators {Ak{^)}. One finds 

Ik,i{^) = 4:Re[{'ip\AkAi\'ip)-{ii\Ak\'ip){i’\Ai\'ip)]. (A7) 

Note that we omitted the explicit dependency of the opera¬ 
tors on the parameters (f. Although the process is unitary the 
QFIM may depend on the parameters Lp, i.e., X = X{Lp). Fur¬ 
ther, in general, we have [A/c((^), Ai{lp)] ^ 0 . 

Appendix B: The QFIM for product probe states 


Next, we define vectors Xn G C^, n = 2,..., 1), with entries 
= {(j)\ak\in)- With this, the QFIM X = 4A^X[^] reduces to 

D D 

I = 4A^ Re[a:„a:|j] = 2A^ [xnX^ + (xnxl)* 

n=2 n=2 

which is a sum of 2{D — 1) rank one matrices. Hence, 
rank[X] < 2{D — 1). 


Appendix C: Unitary multi-parameter estimation: 

Saturating the quantum Cramer-Rao bound 

Next, we prove that the quantum Cramer-Rao bound can 
be saturated in the setting we are considering. In a multi¬ 
parameter estimation setup, in general, the SFDs do not com¬ 
mute. This is the reason why the quantum Cramer-Rao bound 
may not be saturated [ 6 , 7] . As we will see, however, if the 
expectation value of the commutator vanishes, i.e., 

{^pcp\LkLi - LiLkl'ipcp) = 0, (Cl) 

the bound can sill be achieved (see also [28]). One finds 

{'^cp\LkLi /4 

= (■*/’! (Ofe|V’)(V’l + IV’)(V’|oi) X 

X (di\'ip){ij\ + |tA)(t/)|6/) i-ip) (C2) 

= {'>P\Ok\'>P}{‘4’\Oi\^) + {'ip\dk\'ip}{‘>P\Oi\^) 

+ ii’lOkOil'iP) + {i’\ol\'ip){'ip\dj\i^). 


In this section of the appendix, we prove an upper bound on 
the rank of the QFIM for separable probe states. Recall 

I = 4Ar2:[il +4N{N - 1)2:[2] ^ (B1) 

where 

xg = Tr[p[2lafe 0 a;] - Tr[p[ila,] Tr[^[ila,] = 0 (B2) 

for product probe states, i.e., states of the form I'lp) = 
where |0) G C^. Note that for these states = |0)(0| such 
that 

4'i = Re[lV[eWafea,]] - 'IV[^Wafc] IV 

= Re[{(l)\akai\(p)] - {(p\ak\(l)){(l)\ai\(p). (B3) 

Now, let 1 = J2n=i l?n)(?nl where |^i) = |(/)). With this 

D 

4^i = y]Re[((/.|4|C„)(C„|an0)] 

n=2 

+ Re[(^!'|afc|(/))((/)|ai|^!')] - {4)\ak\(t>){(t>\ai\(j)). 

D 

= y^Re[((/.|afc|^„)(^„|a;|<^)]. (B4) 

n=2 


With this. 




= 8 i 




(C3) 


= 8 ilm {'ip\AkAi\'ip) 


where Ok = —lAk and {‘ip\Ak\‘>p) G IR as A*; = For 
Ak = '^n=i expectation value reduces to 


8iy^im[iv[^["’™i4”'®d4'; 

n^m 

= Si^^Im Tr 


-h 8 i^^ Im Tr 




= 8iN Im 


I^Tr , 


(C4) 


since (g) G [R forn 7 ^ m and the last equa¬ 

tion is valid for permutational invariant systems. 

Next, we prove that 


Im 


{'ip\AkAi\'ip) 


= 0 


(C5) 
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is a sufficient condition for the Cramer-Rao bound to be sat¬ 
urated. First, note that each SLD Lk (see Eqn. (A6)) is of 
rank 2 where the non-zero eigenvalues are given by 

= ±2^{i,\Al\i,) - {^l>\AuW (C6) 

with the corresponding eigenvectors 

\c|>^)=auUAk\^l>)+h^U\^), (C7) 


where 


dk 


b 


± 

k 


1 

V2(V’l4lV’) ’ 
(V’lifelV’) ± i\/ 


Hence, the eigenspaces of {L^} are spanned by the d + 1 
vectors 


\^o) = U\^). \^k)=UAkm for k = l,...,d. (C8) 


Secondly, we show that these vectors are linearly independent, 
i.e., the subspace resulting by combining the eigenspaces of 
the SLDs is of dimension d -f 1. To prove this assertion, let 
G G X (d+i) Gramian matrix of the vectors given 

in Eqn. (C8), i.e., Gk,i = (616)- One finds 

/ } ••• ii’lM'iP) \ 

ii’lMi’) ••• ii’lAiAdl'ip) 

^ . . .5 V'^^/ 

XitplAdlip) (iplAdAiltp) ... (iplAdAdlip)/ 

where, of course, the probe state \'ip) is normalised. It remains 
to show that the Gramian matrix has full rank. For this, recall 
that for every Hermitian matrix M that can be partitioned as 

(CIO) 


where A and G are square matrices, it holds that [41] 

M >0-^ A>0andC - B'fA-^B >0, (Cll) 

where M > 0 denotes positive definiteness, i.e., {x\M\x) > 0 
for all \x). Note that S = G — is called the 

Schur complement of block A of M. Now, let A = 1, 
Bk = ('0|A/c|'0),andC/e,z = {i;\AkAi\i;). Obviously, A > 0. 
Further, the Schur complement is given by 

Sk,i = (V’l^fe^lV’) - {'>P\Ak\'iP){'iP\Ai\‘>p), (C12) 

i.e., S = X((^)/4 given that the expectation values of all com¬ 
mutators of the SLDs vanish, i.e., lm[{ 2 p\Aj^Ai\ 2 p)] = 0, see 
Eqns. (A7) and (C5). As we assume that the QFIM has full 
rank (with positive eigenvalues), we have S > 0. Thus, the 
Gramian matrix is positive definite and, hence, has full rank 


such that the set of vectors given in Eqn. (C8) is linearly in¬ 
dependent. Hence, one can find an orthogonal basis of the 
subspace spanned by the eigenvectors of all SLDs by a Gram- 
Schmidt orthogonalisation procedure starting with the vector 
l^o) = U\ 2 p). The d 1 projectors onto these orthogonal 
vectors, together with one element that accounts for the nor¬ 
malisation, form a set of POVMs of cardinality d + 2. As 
one element of this POVM is the projector onto the time- 
evolved probe state, the results of Ref. [8] prove that this set 
of POVMs saturates the quantum Cramer-Rao bound. 


Appendix D: d > 3 

Let us restrict to two-level systems and come back to the set¬ 
ting where the task is to estimate the three components of a 
magnetic field pointing in an arbitrary direction, i.e., the evo¬ 
lution under the Hamiltonian 

3 3 3 

h = flB = J2 f^k^k ^ ^ ^ ^ ^k^k • (f^l) 

k=l k=l k=l 

It is worth mentioning that d = 3 is the maximal number 
of to-be-estimated parameters given the Hamiltonian acts on 
each site independently. Assume that 

d 

h = '^(pkhk (D2) 

k=l 


for a d > 3. We can always decompose each hk in the (nor¬ 
malised) Pauli basis {Pi} with A = dilV2, A = <3-2 /a/2, 
A = ^3 / a/2, and A = 1 / a/2 . One finds 


such that, in fact, [ci] are the independent parameters (and 
the parameters {^k} are determined by the {ck}). Fur¬ 
ther, any contribution that is proportional to the identity can 
be neglected as this would result in an unobservable global 
phase. Hence, estimating these three phases can be interpreted 
as single-particle Hamiltonian tomography at the Heisenberg 
limit. 


4 

Pi = Y,ciPi (D3) 


4 / d 

^ = L] I y^^ipktr[Pihk] 

i=l \k=l 


Appendix E: Single-parameter estimation and multi-parameter 
estimation with commuting generators 

Let us first review the results for single-parameter estima¬ 
tion [4] in the framework discussed in the main text. For this, 
let the single particle Hamiltonian governing the time evolu¬ 
tion be given hy h = iphi^h 2 where the Hermitian operators 
hi and /i 2 do not necessarily commute. Note that this includes 
the estimation of one direction of a magnetic field pointing 
in an arbitrary direction where the remaining directions are 
kept constant, e.g., h = ^ h 2 with ir[cFxh 2 ] = 0 and 
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^ 2 ] 7 ^ 0. As we allow to probe the magnetic field with N 
particles simultaneously, the unitary evolution is given by 

N N 

u = = n (El) 

n=l n=l 

with H = + ^ 2 ) the A^-particle Hamiltonian. 

Hence, A = Yln=i where 

= / da (E2) 

Jo 

such that 


and [Hk^ Hi] = 0 for all /c, / = 1,..., d. For this, we review 
the results obtained in [ 8 ] in the framework discussed in the 
main text. Recall that in [ 8 ] the task is to estimate d phases 
in a d + 1-mode interferometer. Each phase is independently 
imprinted on the probe state in one mode of the interferometer, 
whereas the remaining mode serves as a reference. This is 
done via the generators Hk = where Nk is the number 
operator for mode k. With this 

U{ip) = . (ElO) 

Further, as [TV/., Ni] =0, one finds Ak{cp) = Nk such that the 
quantum Fisher information matrix is given by 




(E3) 


l-k,l 


4Re 




. (Ell) 


Now, for product probe states of the form \ 2 p) = l^n). 

one finds 

N 

1 = 4^ ((,/.„|(aN)2|^„) - (0„|aN|0„)2) (E4) 

n=l 

which reduces to 

I = m {{cp\a^\(l)) - {cp\a\cpf) (E5) 

given that 10^) = |0) for all n = 1,..., A^. The latter is 
maximised by states of the form | 0 ) = + | 0”^^“))/\/2 

where are the eigenstates of d corresponding 

to the minimal Amin(fi) and maximal Aniax(^) eigenvalue. With 
this, 

I = A^(An,ax(a) - An,i„(a))^ (E6) 


The probe state for this QFIM presented in [ 8 ] results from 
the same intuition as the probe state discussed in the main 
text for the magnetic field estimation: is a superposition 

of the states that yield a quantum advantage when estimating 
the parameters individually. While we cannot present a proof 
that this intuition is optimal, it seems a good first guess when 
considering simultaneous multi-parameter estimation. 

Finally, as the generators {Nk} commute, one finds 
(tAliVfeVdV;)* = {i;\NkNi\^P) such that Im [{^P\NkNi\i;)] = 
0 and the Cramer-Rao bound can be saturated. 


Appendix F: Reduced density matrices of the probe state 

In this section of the appendix, we discuss the reduced density 
matrices of the probe state given by 


Allowing for entangled probe states |V^), it is well known that 
the maximal quantum Fisher information is obtained by using 
GHZ-type states [4], i.e., 

+ |</,”“)®^]/^/2. (E7) 


IV;) =V(|$i) + |$2) + |$3)) 

= A4V2 (|$i) + |$2) + |$3)) 

= A4( iv>+)®^+iv.r)®^ 


Note that {| 0 max^(g)JVj 0 min^( 8 )iV| eigcustates of A 

corresponding to its maximal and minimal eigenvalue, i.e., 
{A/^Aniax(^) 5 -A/^Aniin(fi)}- With this, 

2:=A^7Amax(a)-A^„(a))2. (E 8 ) 

Moreover, the Cramer-Rao bound can always be attained 
yielding the quantum advantage of a Heisenberg scaling in 
contrast to the shot noise limit with respect to the precision of 
the parameter cp. Note that for h 2 = 0 and hi = dk, with 
either k = 1, 2, or 3, this reduces to the scenario of estimat¬ 
ing the magnetic field when the direction (here X, Y, or Z) is 
known. 

Next, let us discuss a setting for multi-parameter estimation 
where the generators {Hk} of the unitary time evolution com¬ 
mute, i.e., where 


d 

U = with H{(p) = ^ (fkHk (E9) 

k=l 


where is the eigenvector of the Pauli operator dk cor¬ 
responding to the eigenvalue ±1 for all /c = 1,2,3 and we 
defined \^k) = (10^ )^^)/a/ 2. First, note that the 

normalisation constant is determined via 


1=M‘^ 


6 + 4 



1 -i 
2 


N 


(F2) 



Hence, Ni f/Vh for ^ oc. Next, let us analyse the 
two-body reduced density matrix. First, note that 


+ =tr\2 [\^k){>^>k\] 

= \ + + I) (F3) 

= i (I2 ® I2 + O-fe ® CTfe) 
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for all/c. Moreover, terms liketr \2 [|^/c)(^z|] scale as 1/2^/^ 
such that forN^oo they vanish. Hence, in the limit N —> 
oo, the two-body marginal of the probe state converges to 

k=i 

Finally, let us note that for N = 8 n, n G N, this is exact, i.e.. 

Appendix G: Derivation of the QFIM 

Here, we calculate the QFIM for the probe state given in 
Eqn. (FI) with one- and two-body reduced density matrices 
= l 2 / 2 andp[^] = l 2 ( 8 )l 2 / 4 + 5 ];^_i d/e(g)(j/c/ 12 , respec¬ 
tively. We begin by noting that Tr [d/c] = 0, V/c, since Pauli op¬ 
erators are traceless and, hence, = Tr [d/^d/] /2. With 

this 


such that the entries of X{(p) are given in terms of the entries 
of the operators 

Wk= [da dp d/e (G4) 

Jo 

in the Pauli basis. To find analytic expression of these opera¬ 
tors, recall that with \\n\p = 1 one has 

3 

0 -i^(Efc=i 'rikCTk) — 0 os[d] 1 — i sin[d] ri/ed/e- (G5) 

k=l 

Now, let 


C = V ‘^1 + *^2 + y’i and rik = ^===== (G6) 

V + ^2 + ^3 

for all /c = 1,2,3 (corresponding to X, Y und Z). We find for 
the operators Wk 


= Tr [p^^^d/e (8) di] - Tr [g^^^di] Tr [p^^^d/] 


1 

12 


3 

Tr [{dm ( 8 ) dm){dk ^ di)] 

m=l 



.771=1 


Tr 


Pmdk 


Pmdl 


1 

6 


Tr [dkdi] 


Wi =ai [l + sinc2[C] + (1 - sinc2[^])(7?f - V 2 - vD]/‘^ 
+ 0-2 [1 - sinc^[^]] r]iT ]2 
+ 0-3 [1 -sinc^[^]] rjiTjs, 

where sinc[^] = sin[^]/^. Further, 


as P/e = dkj \/2, k = 1,2,3, together with P 4 = 1 /\/2 is an 
orthonormal basis and the contribution proportional to F 4 for 
the operator d/e is zero. Thus, 




g^^^akai 


= (Gl) 


W2 =(Ti [1 - sinc^[^]] 7]ir]2 

+<72 [1 + sinc2[C] + (1 - sinc2[C])(-7?i + vl - 7?|)] /2 

+(73 [1 - sinc^[C]] r]2T]3, 


and 


Hence, the QFIM is 

4;i = ^ [«,,]. (02, 

Using the definition of the operators {d/e}, see Eqn. ( 6 ), we 
have 


Tr [d/ed/] = [ da dp Tr 
Jo 

Jo 


= da dp Tr 


e‘“'* (7fe ai 




(G3) 


= Tr 


aiWk 


W 3 =CTi [1 - sinc^[^]] ?7 i7?3 
+(72 [1 - sinc^[C]] r/ 2 % 

+ 0-3 [1 + sinc2[^] + (1 - sinc^[C])(-r?^ - vl + vD] /2- 
With this, the QFIM simplifies to 

Ik,i = ^N{N + 2)[{l-smc^[^])rikm + h,isinc^[^]] . (G7) 



















